Reaction cross section described by a black sphere approximation of nuclei 
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We identify a length scale that simultaneously accounts for the observed proton-nucleus reaction 
cross section and diffraction peak in the proton elastic differential cross section. This scale is the 
nuclear radius, a, deduced from proton elastic scattering data of incident energies higher than ~ 800 
MeV, by assuming that the target nucleus is a "black" sphere. The values of a are determined so 
as to reproduce the angle of the first diffraction maximum in the scattering data for stable nuclei. 
We find that the absorption cross section, ira 2 , agrees with the empirical total reaction cross section 
for C, Sn, and Pb to within error bars. This agreement persists in the case of the interaction cross 
section measured for a carbon target. We also find that ^/3/5a systematically deviates from the 
empirically deduced values of the root-mean-square matter radius for nuclei having mass less than 
about 50, while it almost completely agrees with the deduced values for A > 50. This tendency 
suggests a significant change of the nuclear matter distribution from a rectangular one for A < 50, 
which is consistent with the behavior of the empirical charge distribution. 

PACS numbers: 25.60.Dz, 21.10.Gv, 24.10.Ht, 25.40.Cm 



The size of atomic nuclei is considered to be well de- 
duced from empirical data for the proton-nucleus elastic 
differential cross section, da c i/dfl, and the total reaction 
cross section, <tr = cry — a c \, where <jt is the total cross 
section. So far, the analysis that respects both data in 
deducing the nuclear size has not been completed in par- 
ticular for proton incident energies, T p , higher than 800 
MeV. Various approximate theories based on optical po- 
tentials have been proposed to reproduce the elastic scat- 
tering data, while they usually tend to overestimate the 
reaction cross section for 800 MeV < T p < 1000 MeV 
(e.g., Ref. and references therein). 

In Ref. 0, we constructed a method for deducing 
the nuclear size by focusing on the peak angle in the 
proton-nucleus elastic differential cross section measured 
at T p > 800 MeV, where the corresponding optical poten- 
tial is strongly absorptive. In this method, we regard a 
nucleus as a "black" (i.e., purely absorptive) sphere of ra- 
dius a, and determine a in such a way as to reproduce the 
angle of the observed first diffraction peak. If we multi- 
ply a by a ratio between the root-mean-square and 
squared off radii for a rectangular distribution, the result 
for stable nuclei of A > 50 shows an excellent agree- 
ment with the root-mean-square radius, r m , of the mat- 
ter density distribution as determined from conventional 
scattering theories so as to reproduce the overall diffrac- 
tion pattern and analyzing power in the proton elastic 
scattering. 

In this paper, we extend such a previous analysis to 
the case of A < 50, and find out a systematic devia- 
tion between ^/3/5a and r m . We next show that the 
present method is effective at explaining the observed re- 



action cross sections for stable nuclei ranging from light 
to heavy ones. In the black sphere approximation of a 
nucleus, where the geometrical cross section can be de- 
scribed by 7ra 2 , a plays the role of a critical radius inside 
which the reaction with incident protons occurs. We find 
that na 2 is consistent with the measured reaction cross 
section. Consequently, the black sphere picture charac- 
terized by a gives a unified basis for describing the reac- 
tion cross section and the elastic scattering data. This 
simple formula for the reaction cross section, given by 
na 2 , does not include any adjustable parameter. This 
is a feature more advantageous than the fitting formulas 
proposed earlier 0,0 on the basis of the A 1 / 3 law. 

We start with evaluations of a for stable nuclei in- 
cluding light elements by following a line of argument 
of Ref. [2j. The center-of-mass (cm.) scattering an- 
gle for proton elastic scattering is generally given by 
d c .m. — 2sin~ 1 (<2 , /2p) with the momentum transfer, q, 
and the proton incident momentum in the cm. frame, 
p. For the proton diffraction by a circular black disk of 
radius a, we can calculate the value of # c .m. at the first 
peak as a function of a. (Here we define the zeroth peak 
as that whose angle corresponds to # c . m . = 0.) We de- 
termine a in such a way that this value of (9 c m . agrees 
with the first peak angle for the measured diffraction in 
proton- nucleus elastic scattering, 8m- The radius, a, and 
the angle, 8m, are then related by 

2pasin{8 M /2) = 5.1356 (1) 

By setting 

r BS = V3/5a, (2) 
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FIG. 1: (Color) res as a function of mass number, A. For 
nuclei of A > 50, we have plotted the values of r m adopted 
in Ref. Q. For the error bars in res and r m , see Ref. 0- 
For comparison we also plot the results for r m derived in the 
following references: For 12 C, in Refs. H |1 Q (T p = 800 
MeV); for 13 C, in Refs. H Q (T p = 800 MeV]: for 16 0, in 

Refs. dOi ( t p = 1000 MeV ); for 28si > in Ref - II ( t p = 1000 
MeV); for 32 S, in Ref. @ (T p = 1000 MeV); for 39 K, in Ref. 
@ (T p = 1000 MeV); for 40 ' 48 Ca, in Ref. 10] (T p = 796 
MeV), in Ref. (T p = 1000 MeV); for 40 > 42 ' 44 ' 48 Ca, in Ref. 
E3 (T p = 796 MeV), in Refs. H H3 (T p = 1044 MeV); for 
46 ' 48 Ti, in Ref. O (T p = 799.3 MeV); for 48 Ti, in Ref. 
(T p — 1044 MeV). For 4 He, the root-mean-square radius of 
the charge density distribution deduced from electron elastic 
scattering data [l4l Il5| is plotted instead of r m , because no 
r m is available for this case. The crosses ( x ) denote the root- 
mean-square matter radii of the point nucleon distributions, 
and the circles (o) denote those folded with the nucleon form 
factor. The dashed curve represents 0.94A 1//3 fm. The dotted 
line shows A = 50. 



we found that at T p > 800 MeV, res, estimated for 
heavy stable nuclei of A > 50, is within error bars con- 
sistent with the root-mean-square nuclear matter radius, 
r m , deduced from elaborate analyses based on conven- 
tional scattering theory. Thus, expression J5J) works as 
a "radius formula." The factor a/3/5 comes from the 
assumption that the nucleon distribution is rectangular; 
the root-mean-square radius of a rectangular distribution 
is a cutoff radius multiplied by -y/3/5. Here we simply 
extend this estimate of res to lighter stable nuclei. This 
extension of the black sphere analysis is reasonable as 
long as the scattering is close to the limit of the geomet- 
rical optics. This condition is fairly well satisfied at least 
for T p > 800 MeV, since a/A, where A = 2ir/p is the wave 
length of incident proton in the cm. frame, is well above 
unity even for 4 He. As we shall see, the values of res 
are systematically smaller than those of r m for A < 50, 
whereas the values of tto 1 for C, Sn, and Pb agree well 
with the proton-nucleus reaction cross section data for 
T p > 800 MeV. 

The values of res, which are obtained for stable nu- 
clei ranging from He to Pb from the proton elastic scat- 



tering data for T p > 800 MeV, are plotted in Fig. 1, 
together with the values of r m estimated in earlier inves- 
tigations. In collecting the data, we have made access 
to Experimental Nuclear Reaction Data File (EXFOR) 
(lif. In deriving res for 3 He, we have used the data 
for T p = 800 MeV O and for T p = 1040 MeV [3. 
Since the peak position is not clear, we do not include 

" l9l | in this analysis^ For 

_M 

23- We re- 



the data for T p = 1000 MeV 
4 He, we have used the data for T p = 1000 MeV |2 



T p = 1050 MeV [2 



and T p = 1150 MeV 



Ca (T p = 1000 MeV) 
46,48 Ti ( Tp = 799 3 

We do not include 
12 C (Tp = 1000 



mark that, for T p = 890, 991 MeV, no peak has been 
identified since the measured cross sections are limited 
to very forward direction For C and heavier iso- 

topes, we have used the following data: 12 C (T„ = 800 
MeV) 0, E3; 13 C (T p = 800 MeV) H El; 12 C 
(T p = 1000 MeV) JI1l1!C (T p = 1040 MeV) [28|; 14 N 
(T p = 800 MeV) \M, E|; 14 N (T p = 1000 MeV) [13; 
!6o (Tp = 800 MeV) [21 E3; 16 (T p = 1000 MeV) 
[H S_ 20 Ne (T p = 800 MeV) jH^Ne (X = 800 
MeV)]£l; 24M g_2V^ 800 MeV) HI 13; 26 Mg 
(T p = 800 MeV) [Hill; 28 SifT p = 1000 MeV) 
32,34gj Tp = 100Q MeV ) [3(||3cj. 39 K (j, = 10Q0 MeV ^ 

[HEl; 40,42,44,48 Ca ^ = jqq MeV) pj. 40 Ca = 

797.5 MeV) H3 ^Ca = 800 MeV)[3J, E3; 40 Ca 
(Tp = 1000 MeV) |13|4aL4 Ca (T p = 1044 MeV) [H; 
42 ' 44 Ca (Tp = 800 MeV) |43h 42 ' 44 Ca (T p = 1044 MeV) 

ll; 48 Ca (T p = 800 MeV)]4i 1 48 

23; 48 Ca (T p = 1044 Me\0|lli 
MeV) 0; 48 Ti (T p = 1044 MeV)U 
the data for 9 Be (T p = 1000 MeV) 

MeV) 20], nor 16 (T p = 1000 MeV) [20j since the peak 
positions are not clear. 

It is interesting to note that res agrees almost com- 
pletely with r m for A > 50, as shown in Ref. 0, while 
it is smaller than r m by about 0.2 fm for A < 50. In 
order to clarify this deviation, we exhibit the difference, 
^bs — r m, hi Fig. 2. The drastic change in the difference 
around A ~ 50 implies a possible change in the form of 
the real nucleon distribution; the rectangular distribu- 
tion as assumed in the present black sphere model may 
well simulate the real distribution at A > 50, while for 
A < 50 the real distribution is quite different from the 
rectangular one in such a way that the portion of the 
real distribution farther than a is relatively large. This 
feature is suggested by the empirical charge distribution 
deduced from the electron-nucleus elastic scattering [J^ , 
which shows a Gaussian-like form rather than a rectan- 
gular one for light nuclei. 

This feature of the nucleon distribution is expected to 
be reflected by size-sensitive observables for which em- 
pirical data are available for stable nuclei ranging from 
light to heavy ones. Such observables include Is states of 
pionic atoms and isoscalar giant resonance energies; the 
isoscalar part of the pion-nucleus optical potential 
and the inertia associated with the resonances are 
related to the nucleon distribution. 
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FIG. 2: The difference, res — r m , as a function of mass num- 
ber, A. The crosses (x) and the circles (o) are calculated 
from the corresponding values of r m in Fig. 1. The dotted 
line shows res = r m and A — 50. 



We can see from Fig. 1 that res, if having its staggering 
smoothed out, behaves as ~ 0.94A 1 / 3 fm. This suggests 
that res arid hence a provides a measure of the reaction 
cross section, since the empirical values are known to 
behave roughly as oc A 2 / 3 [46J. It is thus interesting to 
compare the black sphere cross section with the empirical 
data. 

We proceed to calculate the proton-nucleus total reac- 
tion cross section from the present black sphere model. 
Our model regards it as the geometrical cross section, 



obs 



na 



(3) 



Here we assume that the incident protons are point par- 
ticles, and that the incident protons, if touching the tar- 
get point nucleon distribution, contribute to the reaction 
cross section by yielding excitations associated with in- 
ternucleon motion. Our model thus predicts vanishing 
cross section for the proton-proton case. This is rea- 
sonable since the proton-nucleon reaction cross section is 
relatively small for T p < 1 GeV HHH. 

By substituting the values of a determined by Eq. 
into Eq. ©, we evaluate obs for stable nuclei. In Fig. 3, 
the results are plotted together with the empirical data 
on a R available for 800 MeV < T p < 1000 MeV 0. 
We can see an excellent agreement between <7bs an d the 
empirical values although the comparison is possible only 
for C, Sn, and Pb. For these nuclei, as shown in Table I, 
<7bs an d <tr agree with each other within the error bars. 
Note that we do not use any adjustable parameter here. 
We thus see the role played by ctbs in predicting <tr, and 
this is useful for nuclides for which elastic scattering data 
are available but no data for an are available. On the 
other hand, if a R is measured, one can deduce a from 
Eq. J2J. In this case, a can be regarded as a "reaction 
radius," inside which the reaction with incident protons 
occurs. In a real nucleus, this radius corresponds to a 
radius at which the mean-free path of incident protons is 
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FIG. 3: (Color) The absorption cross section, <tbs, of a proton 
of T p > 800 MeV by a target nucleus of mass number A. For 
comparison, we plot the empirical data for the proton-nucleus 
total reaction cross section, an (o), which are listed in Ref. |4q| 
for 9 Be, 27 A1, C, Cu, Sn, and Pb. For the latter four elements 
the value of an is the average over the isotopic abundance 
in a target. For these data, we set A as the mass number 
of the most abundant isotope and assign the uncertainty in 
A due to the natural abundance. The dashed curve denotes 
(5/3)tt(0.94A 1/3 ) 2 fm 2 . Inset: cr B s vs a R for C, Sn, and Pb. 
The dotted line represents ctbs = au- 



di the order of the length of the penetration. 

It is natural to generalize the definition of obs given by 
Eq. 10 to the case of the nucleus-nucleus reaction cross 
section. We simply set 



(Tbs = 7r(ai + a 2 ) 2 



(4) 



where ai(a 2 ) is the black sphere radius of a projectile 
(target) nucleus, which we determine from the proton 
elastic differential cross section data for T p > 800 MeV. 

In the case of nucleus-nucleus collisions, rather than 
to measure the reaction cross section ctr, it is more con- 
venient to measure the interaction cross section, aj = 



TABLE I: Values of an and obs f° r C, Sn, and Pb. The values 
of an |46j| are for natural targets, while those of ctbs are for 
12 C, 120 Sn, and 208 Pb, the elastic scattering data for which 
are taken from the references shown in the last column. 



Target 


a R [mb] (T p [MeV]) 


o-BS [mb] (T p [MeV]) Ref. 


C 


209 ± 22 (860) 


217.003 ± 3.45 (800) [6] 
227.788 ± 3.46 (800) [25J 
228.233 ± 9.20 (1000) [2JJ 
237.562 ± 6.69 (1040) [28] 


Sn 


1100 ± 30 (860) 


1117.486 ± 49.6 (800) [49] 


Pb 


1680 ± 40 (860) 


1723.590 ± 88.5 (800) [6] 
1721.680 ± 63.6 (800) [50] 
1606.919 ± 173.0 (800) [49J 
1701.161 ± 94.1 (1000) [30] 
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FIG. 4: (Color online) The absorption cross section <tbs 
for a projectile stable nucleus of A < 30 and a 12 C tar- 
get. We fix r B s = 2.086 ± 0.05 fm for 12 C. For compari- 
son, we plot the interaction cross section, ai (o), measured at 



E/A > 800 MeV for a projectile of 3,4 He, 6 ' 7 Li, 



14,15j^ 16,17,18 r 
12 1 



'Be, 

24,25t\ 



B, 



s O, la F, 2U ' 2i Ne, 2d Na, and 24 ' 2b Mg in- 
cident on a 12 C target [5ll |. The dashed curve represents 
(5/3)tt(2.086 + 0.94A 1/3 ) 2 fm 2 . Inset: ctbs vs <j/ for 4 Hc, 
12,13 C) w Nj i6 0j 20 N6) and 24 Mg The dotted line represents 



0"BS = O-J. 



or — Cincia, where <Ti n eia is the cross section for inelastic 
channels. It is interesting to compare the measured val- 
ues of 07 with the corresponding values of ctbs given by 
Eq. (0J. The results are given in Fig. 4, together with 
the empirical data on 07 measured with a 12 C target for 
stable nuclei of E/A > 800 MeV H^j. 

We find from Fig. 4 that ctbs agrees within error bars 
with the empirical values of 07 except for a few cases. 
This result reassures the role played by the black sphere 
radius as a reaction radius. We also note the tendency 
that 0bs is larger than aj. This is consistent with the 
facts that an > 0/ and that, for a 12 C projectile, 0bs is 
much closer to the empirical value of an |52| than that 

of 0/. 

It is also interesting to evaluate 0bs for highly neutron- 
enriched nuclei once the proton elastic differential cross 
section is measured at high incident energies. If the pro- 
ton reaction cross section 07? or the interaction cross sec- 
tion 07 is measured simultaneously, one could compare 
the result with 0bs- If cbs is significantly smaller than 
0jj (/) , it would suggest that the significant contribution 
to a r (/) comes from the tail region farther than a. This 
could imply the possible presence of a neutron halo |5l| . 
The comparison between 0bs and an (/) provides an im- 
portant check of the widely accepted speculation that one 
can deduce the halo structure from measurements of the 
interaction cross section 

In summary we have extended our previous analysis 
of the nuclear size for A > 50 2] , which is based on the 
black sphere assumption of a nucleus, to lighter stable 
nuclei. The black sphere radius, a, has been determined 



so as to reproduce the angle of the first diffraction maxi- 
mum in the proton elastic scattering data. We have thus 
finalized a systematic analysis of the existing data for 
proton elastic scattering off stable nuclei ranging from 
He to Pb at proton incident energies above ~ 800 MeV. 
We have two significant results. First, the values of rss 
obtained from Eq. J5J) are systematically larger than the 
root-mean-square radius r m deduced from elaborate scat- 
tering theory for A < 50, although they agree quite well 
with each other for A > 50. This suggests a significant 
deviation of the nucleon distribution from the rectangu- 
lar one for A < 50. Second, the absorption cross section, 
7ra 2 , is consistent with the empirical total reaction cross 
section for C, Sn, and Pb. This consistency persists in 
the case of the interaction cross section measured for a 
carbon target. We thus see the dual role of a as a black 
sphere radius and as a reaction radius. 

Elastic scattering and total reaction cross section data 
for heavy unstable nuclei are expected to be provided by 
radioactive ion beam facilities, such as GSI, Darmstadt, 
and the Radioactive Ion Beam Factory in RIKEN. Ex- 
periments at RIKEN are planned for a beam energy of a 
few hundreds MeV per nucleon. It is thus important to 
study the energy dependence of the black sphere radius, 
a. As a first trial, we have studied the case of 208 Pb, for 
which not only the de Broglie wavelength of the incident 
proton but also the mean-free path of the proton in the 
nuclear interior is sufficiently short compared to the nu- 
clear radius, and confirmed that the resultant obs agrees 
with the measured value of an within the error bars for 
the proton incident energy down to a few hundreds MeV. 
The energy dependence is found to be similar to that of 
the nucleon-nucleon total cross section. This is consistent 
with the fact that the black sphere radius corresponds to 
the reaction radius. We remark that the energy depen- 
dence of a may well modify the relation between a and 

Towards future possible application of the black sphere 
model to neutron-rich unstable nuclei, we here give a cou- 
ple of examples, 6 ' 8 He and 11 Li, for which empirical infor- 
mation on the proton-nucleus total reaction cross section 
is available at relatively high energy. Hereafter we will 
estimate the black sphere radius a from an . — na 2 in the 
absence of the data for the first peak in the elastic diffrac- 
tion pattern. From the empirical data, an = 161.3 ± 3.7 
mb for 6 He (E/A = 721 MeV) and a R = 197.8 ± 3.5 mb 
for 8 He (E/A = 678 MeV), measured by Neumaier et al. 
0, we estimate a = 2.27±0.03 fm and r BS = 1.76±0.03 
fmfor 6 He, and a = 2.51 ±0.02 fm and r BS = 1.94 ±0.02 
fm for 8 He. The values of r m deduced from the elastic 
scattering data are r m = 2.45 ± 0.10 fm for 6 He and 
r m = 2.53 ± 0.08 fm for 8 He. We thus see a difference 
between r m and tbs of order 0.7 fm, which is significantly 
larger than that for light stable nuclei (see Fig. 2.) For 
11 Li, we take note of the data for the interaction cross 
section obtained for a proton target: aj = 276 ± 8 mb 
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at 800 MeV/A 57]. If we assume ctr = 07, we obtain 
a = 2.96 ± 0.04 fm and r BS = 2.30 ± 0.04 fm. Since 
the typical value of r m amounts to about 3.1 fm |5l| . 
we obtain r m — res of about 0.8 fm, which is similar to 
the case of neutron-rich He isotopes. This leads to an 
interesting implication that such a difference is typical of 
nuclei having a neutron halo. In order to justify this im- 
plication, however, it is inevitable to confirm, by future 
experiments, the assumption that even for neutron-rich 
unstable nuclei, (Jr agrees with the black sphere cross 
section determined from the first peak angle in proton- 
nucleus elastic scattering differential cross section. 
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